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SUMMARY 


The conventional six-engine reaction control jet relay attitude 
control law with deadband is shown to be a good linear approximation 
to a weighted time-fuel optimal control law# Techniques for evaluating 
the value of the relative weighting between time and fuel for a 
particular relay control law is studied along with techniques to 
interrelate other parameters for the two control laws# 

Vehicle attitude control laws employing control moment gyros 
are then investigated. Steering laws obtained from the expression 
for the reaction torque of the gyro configuration are compared to a 
total optimal attitude control law that is derived from optimal linear 
regulator theory# This total optimal attitude control law has com- 
putational disadvantages in the solving of the matrix Riccati equation. 
Several computational algorithms for solving the matrix Riccati 
equation are investigated with respect to accuracy, computational 
storage requirements, and computational speed# 
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I. INTRODUCTION 


Some aspects of optimal attitude control systems are examined in 
this study. The conventional six-engine reaction control jet relay 
attitude control law with deadband is compared to an optimal weighted 
time-fuel attitude control law. The relay attitude control is shown 
to be a good linear approximation to the weighted time-fuel optimal 
control law, then a procedure is developed to determine the optimal 
control law weighting of time and fuel that corresponds to a given 
relay attitude control law with deadband. This time-fuel weighting 
value is determined by selecting the value that minimizes the mean 
square error between the switching curves of the two control laws. 

Next, a procedure is developed to determine the six-engine relay atti- 
tude control law with deadband that most closely corresponds to a 
given optimal weighted time-fuel attitude control law. 

Vehicle attitude control laws employing control moment gyros 
(CMGs) are then investigated. The reaction torque of the gyro con- 
figuration on the vehicle is determined in terms of the gimbal angles , 
the gimbal rates, and the vehicle rates. This expression cannot be 
solved for the gimbal rates directly, so different procedures for 
obtaining the gimbal rates are examined. Basic approaches to obtain- 
ing these steering laws given the commanded torque are compared with an 
optimal control moment gyro attitude control law that is obtained 
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from optimal linear regulator theory. This system is reduced to the 
linear regulator by frequently linearizing the system as the CMG 
gimbal angles travel some prescribed amount from their previous 
nominal value. This total optimal CMG attitude control law has some 
computational disadvantages in attempting to solve the matrix Riccati 
equation on a flight control computer. A sensitivity analysis is 
performed on the gain matrix to determine if any gains can be considered 
constant. Several computational algorithms for solving the Riccati 
equation are investigated with respect to accuracy, computational 
storage requirements, and computational speed. 



II. RELATION OF A REACTION CONTROL JET CONTROL LAW 
TO AN OPTIMAL WEIGHTED TIME- FUEL CONTROL LAW 

The attitude control law for the conventional six engine reaction 
control jet system has the phase plane portrait shown in Figure II-l. 
Although it is suspected that this control law was originally developed 
empirically, it will be shown that this law closely approximates one 
determined by utilizing optimal control procedures. 

The analysis will be simplified by assuming that the system 
dynamics can be represented as 1/s • It will also be assumed that the 
position, <fp, can be obtained by integrating the body rate, <f> » directly 
and that the position and rate signals can be measured. With these 
simplifications, the system to be controlled is depicted in Figure II-2. 
The state equations for this single-axis system are 



(II-l) 


For a reaction control jet control system, the elements to be 
considered in a performance measure are time, fuel, and final state. 
The performance measure will not penalize the final state, but rather 
the final state will be constrained to lie within the deadzone limits 
(i.e., | <t>(tf ) | ^(Jb) * Thus, the performance measure weighting time and 
fuel is 
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Figure II-l* Phase plane portrait of six-engine 

reaction control jet attitude control law. 




Figure II-2. Open loop simplified system to be 
controlled. 
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J(u) = 



* + |u<t)| 


> dc 


(II“2) 


where \ is the relative weighting of time and fuel, and u(t) is the 
control law. Since the final position is constrained to lie within 
the deadzone limits, this problem can be solved for $(tf) = 0 and the 
results shifted to the limits and <|>(t^) = for the 

solution to the deadzone problem. 

For this problem, the Hamiltonian is 

H ~ \ + | u | +p^*2 + P2 U (H-3) 

The necessary conditions for optimality are 

dH_ m St L m & 2 (II-4) 

9pi 

• 

M m st 2 m a (n-5) 

3p2 

« 

- ML = Pi = 0 (n-6) 

-M_-P 2 = -h (II-7) 

9x 2 ” 

M + P 2 a <_ juj + P 2 u (n-8) 

where the hat above the states, costates, and control indicates the 
optimal trajectories and the optimal control history, Pontryagin's 
minimum principle. Equation (II-8) , reveals that the form of the optimal 
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control is 


a 


/+1 for p 2 < “1 

0 for -1 < f> 2 < 1 
/ * 

-1 for 02 >1 
\ singular for p 2 “ ±1 


(IX-9) 


Integrating Equations (II-6) and (II-7) yields the following solutions 
to the costate equations 


h - c i 

02 = “ c l t + C 2 


(11-10) 


where and C 2 are constants of integration. The solution for £2 i* 1- 
dicates two things about the optimal control since fl is expressed as a 
function of the costate trajectory, £2* First, singular control exists 
only for the special case with C^=0 and C2=±l. If is any value 
other than zero, £2 cannot remain at ±1 for any finite period of time. 

For C^=0, 02 =C 2 indicatin S the condition that C2=±l for the existence 
of singular control. Secondly, for C^O, p 2 changes sign no more than 
once (at 1=02/0^) and passes through P2 =- 1 f>2 = +l ao more than once, 

resulting in two control switchings at most. This indicates that the 
form of the optimal control law in its most general form will be 
{1,0,-lJ or {-1,0,1}. Depending on the initial conditions, the first 
elements in these control strategies might be unnecessary. These special 
case control strategies are {0,-1}, {0,1}, {-1}, and {1}. 
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Integrating Equation (II-5) with d=±l yields the solution for 5^ as 


= ±t + C 3 


( 11 - 11 ) 


where C 3 is a constant of integration. 
(II-4) yields 

% = ±ht 2 + c 3 t + c 4 


Then integrating Equation 


(II-12) 


where C 4 is a constant of integration. Solving Equation (11-11) for t 
and substituting this result into the expression for given by 
Equation (11-12) gives the trajectories in the phase plane for u=±l as 

*1 = f C 5 (11-13) 

where C 5 - C 4 - J* ,c| for u = +1 

and = C 4 + for u = -1 

Integrating Equation (II-5) when d=0 reveals that $2 is simply a constant 
for any such period in the optimal control history. Then, integrating 
Equation (II-4) yields 


^1 “ ^t + ^2 


(11-14) 


where k^ and k -2 are the constants of integration. Therefore, when ti=0, 
#2 remains constant and increases or decreases with time depending on 
the value of &2 a t the time d is switched to zero. 
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The parabolas defined by Equation (11-13) for C^=0 are the min- 
imum time switching curves that will bring the system into the origin of 
the phase plane. The remaining step is to determine the switching lines 
that define the switch from u=-l to u a 0 and the switch from h=+l to h=0. 
Following the procedure employed by Kirk Jl] , t Q denotes the time when 
the optimal control switches from +1 to 0 and t^ denotes the time when 
the optimal control switches from 0 to -1. Figure 1 1-3 shows some 
candidate trajectories for a given initial condition. Since t Q occurs 
somewhere on the segment C-K, points D, F, and H are candidate points 
on the switching curve that switches 0 from +1 to 0. The points E, G, 
and I are corresponding candidate switching points for the switch from 
0 to -1 since t^ must occur on the segment K-0, Equation (11-13) relates 
and &2 on K-0, so that 

%(t) = -*5*2 (t), (11-15) 

which leads to 

%(t x ) = -^( t ± ) . . (11-16) 

Integrating Equation (II-5) with fl=0, for the switch from d=+l to ti=0 
occurring at t=t 0 , yields 

* 2 (t) = k 3 = $(t o ) (11-17) 

Then, integrating Equation (II-4) gives 





10 


CH . ( t 1 

%(t)dt - S 2 (t o )dt 

) t Jr 


(11-18) 


Therefore, performing the integration of (11-18) yields 

S l (t l> = % (t o ) + *2 (t o )It: l " t o ] (H-19) 

As was earlier determined using Pontryagin’s minimum principle, f>2 =± ^ 
for control switchings. Using this result in Equations (11-10) yields 

^ 2 (t o ) = "Vo + C 2 “ "I (H-20) 


W = - c l c l + C 2 - 1 


(H-21) 


Since £ 2 (t) is a constant for 11=0, then £ 2 (t 0 ) = it 2 ( t i)* Since the 
Hamiltonian is required to be identically zero along the optimal 
trajectory, at t=t Q 

X + |<l(t 0 )|+ C 1 ± 2 (t (> ) + 0 2 (t o )u(t o ) = 0 (11-22) 

but |fl(t o ) j + ^ 2 (t 0 )d(t 0 ) ° 0 (11-23) 

since = Therefore, 

X + Cj* 2 (t 0 ) = 0 (11-24) 

Solving Equation (11-24) for gives 
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C 1 = - * (11-25) 

*2(t 0 ) 

Subtracting Equation (11-21) from Equation (11-20) yields 

Itl - t Q ] = -2/q , (11-26) 

which, using Equation (11-25), becomes 

Iti - t Q J = 2x 9 (t n ) (II-27) 

\ 

i 

Equating the expressions for ^(t^) given by Equations (11-16) and 
(11-19) gives 

-*5*2 < t x ) = -*sA|(t 0 ) = it 1 (t 0 ) + (t 0 ) rt x - t Q ] (H-28) 

Using Equation (11-27), Equation (11-28) becomes 

-^(t 0 ) = %(t Q ) + 2Ag(t r ) (11-29) 

X 

which can be solved for ^(t Q ) to give the switching line in the phase 
plane as 

*lCt ) = „ *iCt Q ) (11-30) 

2X 

In an analogous manner the switching curve for 0=-l to 0=0 is determined 


to be 
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= l±L *i(tJ (11-31* 

2X 

Figure II-4 shows several of these optimal switching curves for 
various values of X . These results confirm the intuitive feeling that 
as X approaches infinity the weighted time-fuel law approaches the 
minimum time control law. Figure II-5 shows some weighted time- 
fuel optimal trajectories for three sets of initial conditions with X=l. 

Now, to analyze the original deadzone problem, these results will 
be shifted to <Ktf) - and 4 >(tf) = - 4 ^ by a change of variables. 

The phase plane portrait that is obtained by this shift is shown in 
Figure II-6. However, since the control law only maintains within 
a deadband, then the inner switching curves are meaningless. The 
internal X switching curves, B-F and C-G, will be omitted because 
another of the switching curves has preceded these curves in performing 
its originally intended switching function. The internal minimum time 
optimal switching curves, A-B-C and F-G-H, could possibly be used 
along with the external switching curves, I-C-D and E-F-J. They are 
omitted because any overshoot resulting from applying u=±l along the 
internal minimum time switching curves places the system out of the 
deadzone. Their omission also more clearly indicates that the optimal 
control is 0=0 in the shaded region of Figure I 1-6. Since this 
contains part of the deadzone region, control effort is unnecessary 
because the control objectives have been met. Since the internal 
minimum time switching lines, B-C and F-G, are special cases of the 
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fuel switching curves for X =.1,1,10. 
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Figure 11-5. Optimal weighted time-fuel trajectories for intial 
conditions A,B,C with X=l. 
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Figure II-6. Shifted weighted time-fuel optimal switching curves to 

give the optimal de adz one weighted time-fuel control law. 
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region just mentioned and 6=0 on both sides of these curves, they will 
be omitted. The resulting optimal attitude control law composed of 
switching curves I-C-D and E-F-J in Figure II-6 is shown in Figure II-7 
along with the presently employed reaction control jet attitude control 
law superimposed for comparison. 

The control system that gives the presently employed reaction 

control jet control law is shown in Figure II-8. Before determining 

an expression for a value of A that gives the optimal weighted time- 

fuel control law that corresponds to a particular rate ledge and rate 

limit control law, some of the terms shown in Figures II-7 and 

II-8 will be defined. <f> R is the limiting value for $ in the (j> , 

feedback loop in Figure II-8. Figure II-7 reveals that the sloping 

deadzone curves reach their limiting values for 4>=±4> . The maximum 

R 

rate of the system in the 0=0 region is while the minimum rate 

in the 0=0 region for is <^ e< jg e * <t> c is the commanded vehicle 

attitude, and is the width of the deadzone about <f> c . Finally, 
is the rate feedback gain. 

The rate ledge and rate limit can be related to the states which, 
in turn, are related to A in the phase plane. These expressions may 
be solved for the value of A that gives a weighted time-fuel curve 
which intersects the rate ledge at <J> = 4 >q 38 shown in Figure II-7. 

For the equations for the switching curves can be determined 


using Figure 11-8 to be 
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Figure II-7. Phase plane of the reaction control jet control law and 
the corresponding optimal control law. 



Figure II-8. Single-axis reaction control jet control system. 
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A]^ + <j> - 


Aj_$ + 4 1 = — 


(11-32) 


Equation (11-13) defines the corresponding optimal switching curves as 


hi 2 + tf> - <^dB 
~hi 2 + $ = -<t>dB 


(11-33) 


Substituting for cj> in Equations (11-32) and (11-33), solving these 
equations for <f>, and equating the expressions for 4 gives the expression 
for to be 


A 1 * (11-34) 

The rate limit and rate ledge are obtained from Equations (11-32) by 
solving for <}> with These expressions for the rate limit and 
rate ledge control law are 


♦lim ~-J 2 (*R+ *dB>' (H-35) 

Pledge - (*R ' *<a\f 2 " <II - 36) 

V 'J'R + 'J’dB 


The expression for the optimal switching curve corresponding to Equa- 
tion (11-36) is given by Equation (11-30). Solving Equation (11-30) for 
and including the shift from the origin to -ij^B gives 

+ 2 * - rh ( * + *dB> 


(XI-38) 
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Squaring Equation (11-36) and equating the expressions for yields 


2 ^R ~ ^dB )2 
+ ^dB 
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1+4 


O + 


^dB^ 


(11-38) 


Equation (11-38) may be solved for either 1 or corresponding to the 
intersection of the rate ledge and the optimal switching curve. Solving 
for $=<f )0 yields 


*0 - - X + 4 (tR ~ fdB ? 2 -‘f’dB 
* ^R + ^dB 


(11-39) 


while solving for 1 gives 


~ 4< ^r ~ W 2 

(4>r “ 4>dB) 2 + (4> + + ^dB^ 


(11-40) 


The optimal weighted time-fuel control law that most closely 
corresponds to a particular reaction control jet (RCJ) control system 
may be determined by minimizing the mean square rate errors between 
the two switching curves. The value of 1 determined minimizes the 
mean square error for a given value of <)> , the largest expected 

attitude error. The square error, f(\,<t>), is defined as 

f »4>) ^ ^optimal ” ^rate ledge^ (11-41) 

Solving Equation (11-37) gives the egression for as 

^optimal “ ~ — — — ^ + ^dB^ ^ (11-42) 
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Pledge *- s obtained from Equations (11-32) and (11-36) as 




rate ledge 


-C* + * dB ) 


[ ( *R ~ W I 


+ ^dB 


> + *d b/ 


" ^dB 


(11-43) 


1/2 » "♦max < * < “ *R 


The value of X that minimizes the mean square error is obtained by 
solving 


3 

- 1 i 

"^max 

f(X ,<f>)d(f> 

ax 

<h - 

r max y dB 





(11-44) 


The integral in Equation (11-44) must be divided into two regions since 
♦ 

<J>rate ledge * s defined by two different expressions over the interval 
from -if 1 ^ to “$ max » Substituting Equations (11-42) and (11-43) into 
Equation (11-41) and performing the squaring operation yields 


where. 


ffr ><j>) - 



’ "♦ R < * < ' ^dB 


* “+max < * < " *R 


,*) = _ 2A (4> ± j> d p ) 4(<f> + 

X + 4 


. U + W- 

(V + 4) (<J> R + <J> dB ) 


(11-45) 


1/2 


J 


2£U_£d&i 2 

^R + ^dB 


(11-46) 
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311 d » 40 = — 2A (j> ^ 4>^g ? 

A + 4 

? C^R — idftl 2 
*R + ^dB 


Using (11-45) , Equation (11-44) now becomes 


- 4 C*r - <j> dB ) ± jj ± 4>an) 


1/2 


Oi + 4) C$ R + 4> dB ) 


(11-47) 



7 ‘ *R 

^max 


r f i ex , <t>)d$ + 


ax 

iL 

), 

“*R J 


- 0 


(11-48) 


where Equation (11-44) has been multiplied by the constant (^ max ~ 4*dB^ 
Evaluation of equations in the form of Equation (11-44) is aided 
with the use of Leibniz rule, which is 


JL_ 

at 


b (t) 


a(t) 


g(t,u)du = 


b (t) 


a(t) 


9g(t»u) du + g(t,a(t)) db(t) 
3t dt 


g(t,b(t)) da(t) 
dt 


(11-49) 


Applying Leibniz rule. Equation (11-48) becomes 


/--*R , ^ , 

C ^max 


\ afjQ^ » 4 *) + 

\ 3f oft ,*> d<j> = 0 

( 11 - 50 ) 

) 3 ( 

J 


-l ^dB 

_ ^R 



The partial derivatives required in Equation (11-50) are obtained by 
taking the partial derivatives with respect to X of Equations (11-46) 
and (11-47) . These partial derivatives are 



, <j>) 


& 


3f 2^ »$) 

3V 
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8 (^> + 4>dfi) , 8 

[-(* + <f> dB .)] 3/2 

(X + 4) 2 

[*' C* E + * dB >] 1/2 .c + 4> 372 

8 (4> + ^>dB) _ 8 ^R " ^dB^[“^ + ^dB )] 1 ^ 2 

a + 4) 2 

! u *r + ♦dB ) ] 1/2 + 4)3/2 


Substituting Equations (11-51) into Equation (11-50) and evaluating 
the indicated integrals yields 


where 


"1 *2 3 

* 1 / 2 < M - 4> 3/2 (\ + 4) 2 X 1 / 2 (* + 4) 3/2 

K _ 8 f <2 /S)C» r - W 5/2 ' 

(4 >r + < J > dB )!/2 


(11-52) 


(II-53a) 


*2 4 ^max “ ^dB^ 2 

K 3 - 8L(2/3K^ - j)Jf (W - +dfi) 3/2 - <*R - *dB> 5/2 ] 
(♦R + *dB) 1/2 


(II-53b) 


(II-53c) 


Co shining like terms in Equations (II-53a) and (ll-53c) , then substituting 
the result in Equation (11-52) gives 


K/ 


1 1 ^ 2 ft + 4) 3/2 


K 2 


& + 4)' 


= 0 


(11-54) 


where K 4 = 0-6/3) 0fr R - <{> dB )(j> maJt - - (32/15) (<fr R - (ft^) 5 / 2 

(♦r + *dB> 1/2 

(11-55) 
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Multiplying Equation (11-54) by X (A + 4)^ and transposing terms so that 
only terms involving A are on one side of the equation yields. 


I l_f /2 = K 4 (11-56) 

A + 4 ) K 2 


Squaring Equation (11-56) and solving for A gives the expression for 
the A that minimizes the mean square rate error between the two switching 
curves as 


4Kf} 

k£ - K$ 


(H-57) 


As an example, suppose a RCJ control system has = 2, <(ig = 4, 
and (J’max = *0. Figure II-7 shows the relationship of these terms to 
the control law. To find the weighted time-fuel optimal control most 
closely corresponding to the RCJ system, the constants K 2 and are 
first computed using Equations (11-53) and (11-55), respectively. 

For the example, these constants are calculated to be 

K 2 = 256 

K 4 = 93 

Substituting these values into Equation (11-57) yields A as 
A = .61 

This system places more emphasis on conserving fuel since the control 
law penalizes fuel consumption 1.7 times as heavily as it does time. 
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However, since fuel is limited to the fixed amount the vehicle can 
carry, a RCJ system will usually place more emphasis on fuel conservation. 

The previous results were obtained from the standpoint of 
beginning with the conventional RCJ control system with the parameters 
4*dB* (j> R , and given. The optimal weighted time-fuel control system 

that most closely corresponds to the RCJ control system is obtained by 
solving Equation (11-57) for the X that minimizes the mean square 
error between the two rate ledge switching curves. Now the problem 
will be approached in the reverse direction. The value of X in the 
weighted time-fuel performance measure along with the desired and 
<J> max for the RCJ control system are given. The problem is to find the 
RCJ system most closely corresponding to this optimal control system. 

An expression for can be obtained in a manner analogous to the 
technique used to obtain an expression for X . For this case, Equation 
(11-44) becomes 


9 


1 

^max " *dB 



$)d<}> 


(11-58) 


where f(4> R » <J>) is the same expression as that defined for f(k , <jj) in 
Equation (11-41) with the roles of X and as known and unknown vari- 
ables interchanged. Applying the same procedure employed to determine 
an expression for \ results in a sixth order expression in This 

expression can be solved with a digital computer using root-solving 
algorithms. However, a more appealing approach is presented here. 
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The initial step is to obtain a second order expression in <J> R in 
terms of 4>g using Equation (11-39) . Rearranging Equation (11-39) and 
performing the squaring operation on (<}> R - W y ields 

i R 2 + a 1 | R +a 0 *° (11-59) 

where a^^ = -2<j> dfi + — (<t> 0 + 

A .+ 4 

a 0 = ^dB 2 + “ ^ — ( *0 + W*dB 
X + 4 

The quadratic formula is used to determine <J> R for a given i^q as 

<f>R = - a l ±(a ? ~ 4a 0 )1/2 (11-60) 

2 

Equation (11-60) is used to obtain an initial guess of the value of $£• 
This value of 4> R is obtained by initially guessing the value of <j>g to 
be midway between (f)^ and <t> max since this would be a reasonable guess 
for the value that minimizes the mean square error. Therefore, 

4>q = -* dB - W 'Adi (11-61) 

2 

for the initial guess of <)>q. 

In the expression for A given by Equation (11-57) , only is a 
function of <j>£. For the given value of X , the desired value of K^, K^, 
is obtained by solving (11-56) to give 
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K,. = K ./ X | 1/2 (11-62) 

4 1^1 

A Taylor series expansion of is obtained from Equation (11-55) in 

order to systematically improve the initial guess at <f> R . Using only 
the linear terms in the Taylor series. Equation (11-55) becomes 


K 4 (* r ) = K^rq) + VUj^) 6 * r (11-63) 

where 6* R = <fr R - ^ 

^'(♦jo ) = (l6/3m max - - (32/15) (C - D) (II-64a) 

B = (^RO + ^dB^ 1 ^ 2 ~ >5 ^R + ^dB^ ^^RO ~ ‘W (II-64b) 

(<*R0 + ^dB^ 

C = 2.5(4>rq - $ dB ) 3/2 (4> RO + W 1/2 (II-64c) 

( $R0 + W 

D = ,5 ^*R0 + ^dB^ l/2(<t> R0 ~ ^dB )5/2 (II-64d) 

(4>R0 + ^dB^ 


<f> R g is the present guess for <f> R and K 4 ((j> R Q) is given by Equation (11-55) . 
Solving (11-63) for 6<f>R used to update iJ> R yields 

6<f> R = K 4 ^r) ~ K4(*R0> = K 4d ~ K 4^RC)) (11-65) 

K 4^^R0^ K 4'0f>R(p 
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The improved guess for <{>r is now given as 


*R = *R0 + 6 *R (11-66) 

As an example of applying this procedure, consider the earlier 
example where \ was computed to be .61 with = 2 and <l> maK = 10. 

The initial value of 4 >q is computed from (IT-61) as 


*0 = " 6 

Using this value of 4 >q , <j> R is determined from (11-60) to be 
<J> R = 3.74 

K^OJiro) is computed from (11-55) as 

K. (3 . 7 4 ) - 83.5 
4 


is determined from (11-62) to be 


K 4d * 53 


Next, solving (11-65) for the update, 5<j) R , gives 


<$<j>p = 93 - 83.5 _ .252 

37.8 

The new value of <j> R is given by (11-66) as 
$ R = 3.74 + .252 = 3.992 
Computing the new value of ^(^^q) gives 


K 4 (3.99) *=92.9 
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For this case only one update was required to obtain an accurate value 
for This value of ^ can then be used to determine using 

Equation (11-34), the rate limit using Equation (11-35), and the rate 
ledge using Equation (11-36). These values are computed to be 

A 1 = 1.73 
♦lim = 3 - 46 
Pledge = 1 * 16 

These values along with the given ^g determine the RCJ control system 
that most closely corresponds to the given optimal weighted time-fuel 
control system. The results obtained in these examples illustrate 
the uniqueness of these procedures. That is, if the procedure is 
applied to find an optimal weighted time-fuel control system correspond- 
ing to a given RCJ system and then the reverse procedure applied to 
this optimal time-fuel system, the RCJ system obtained will be the 
original one. Figure II-9 shows the optimal control system and its 
corresponding RCJ control system for the two examples given. 

In summary, under the assumption of second order dynamics, the 
conventional RCJ control system is shown to be a good linear approxima- 
tion to an optimal weighted time-fuel control system. Furthermore, for 
4>dfl> Lav i and ‘f'R given, an expression is derived that gives the value 
of X weighting time and fuel. This A is chosen such that the mean 
square error between the switching curves of the RCJ system and the 
optimal control system corresponding to this value of A is minimized. 
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Also, an iterative procedure is presented to determine for the RCJ 
control system that corresponds to an optimal weighted time-fuel 
control system with X given. This procedure also requires knowledge 
of the desired and for the RCJ control system. The procedure 

produces unique results in that applying the procedure followed by 
its reverse procedure yields the original system. 



^max = 

Figure II-9. Phase Plane for Example System. 






III. SOME BASIC IDEAS RELATED TO THE USE OF CONTROL 
MOMENT GYROS FOR ATTITUDE CONTROL 


A reaction control jet (RCJ) mass expulsion system can provide 
control within a band of a desired operating state. Control moment 
gyros (CMGs) , based on the principle of momentum exchange, can provide 
the fine-pointing capability required for many space missions. Morine 
and O'Connor describe the CMG and its relative merits in {2]« 

The principle of momentum exchange is a consequence of Newton's 
second law of angular motion which states that the total external 
moment, M exfc , acting on a system, is proportional to the time rate of 
change of angular momentum with respect to inertial space. This can 
be written as 




Inertial Space 


(HI-1) 


Considering the system to be composed of a vehicle and a controller. 
Equation (III-l) can be integrated to yield 

J”ext dt ’ "ext * “control + “vehicle ’ “system®*) < m - 2 > 
Equation (III-2) illustrates that the controller momentum can be used to 

i 

balance external torques as well as change the spacecraft attitude by 
varying the vehicle angular momentum. 
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Since electrical energy is used as the prime source of power, 
the concept of momentum exchange using CMGs is desirable because it 
provides continuous vehicle control and a recoverable energy source. 

Cyclic disturbance torques can also be handled on a continuous basis 
over a long period of time. However, constant external disturbances 
applied to the vehicle will eventually cause the controller to reach its 
maximum capacity, and thus saturate the CMG. This saturation will re- 
quire the expulsion of propellants to remove some momentum from the 
saturated CMG. This propellant expulsion task may be accomplished at 
a convenient time during the mission. 

The CMG is essentially a gimballed wheel rotating at a constant 
speed which provides a constant angular momentum magnitude capable of 
variable orientation relative to the spacecraft. Exchange of momentum 
between the vehicle and the controller is effected by causing a change 
in direction of the constant momentum magnitude of the CMG. 

Figure III-l shows schematically a two degree of freedom CMG. 

It consists of a wheel that rotates at a constant speed. This wheel is 
held in a housing which is called the inner gimbal. The inner gimbal 
is coupled to the outer gimbal through the Cl) pivot which is perpendicular 
to the wheel spin vector as shown in Figure III-l. The outer gimbal is 
held to the base through the (3) pivot which is perpendicular to the 
(1) pivot. The (1) and (3) pivots are driven by geared D. C. torque 

motors. 

CMGs offer several advantages over other momentum exchange 
devices such as reaction wheels. Since the CMG operates at a constant 




Figure III-l. The "j M control moment gyro schematic. 
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speed for which the efficiency can be optimized, its efficiency is 
greatly improved for most operating momentum values. Larger maximum 
moment control can be obtained will relatively small turning rates of 
the large constant momentum. Another CMG advantage is better bandwidth 
over the reaction wheel which has a limited bandwidth due to the physi- 
cal characteristics of the wheel. A final advantage of the CMG is that 
it has a more linear approach to saturation than do reaction wheels. 

These advantages coupled with the fact that the CMG provides continuous 
control and a recoverable energy source makes the CMG an attractive 
device for long duration fine-pointing attitude control missions. 

A. CMG Equations of Motion 

For any control System using CMGs, it is first necessary to know 
the relationship between the reaction torque of the CMG configuration 
and the gimbal angles and gimbal rates of the gyros. Throughout this 
analysis , the CMG configuration ass ume d will be the conventional S1XPAC 
configuration as shown in Figure III-2. 

The relationship between the gyro gimbal angles and the reaction 
torque will be obtained by first arriving at an expression for the 
angular momentum of the configuration, and then applying Equation (III-l) 
to determine the torque. 

The initial step in determining an expression for the angular 
momentum is to define three gyro spaces [ 2 ] . Referring again to Figure 
HI-1, the space or inner gimbal is described by the coordinate 

system l 1A(j) > l 2 A (j) ’ andl 3A(j)‘ The A (j) 


space coordinate system 



OUTER ACTUATOR 
ASSEMBLY 



1C(1) 


INNER ACTUATOI 
ASSEMBLY ^ 



l 3C(3) 



CMG 3 


INNER ACTUATOR ASSEMBLY 

OUTER ACTUATOR 
/ASSEMBLY 


CMG 2 


1C (3) 


CMG 1 

OUTER 
/ACTUATOR 
/ ASSEMBLY 



3C(1) 


INNER 

ACTUATOR ASSEMBLY 


Figure III-2. The Conventional CMG SIXPAC configuration 
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is defined such that the l2A(j) vector is along the spin reference 
axis, the vector is along the Cl) or inner pivot, and the l^ACj) 


vector completes the orthogonal coordinate system. The ”C 


ftp " 


0 ) 


space or outer gimbal space is described by the coordinate system 


1 lC(j ) * 1 2CQ)* 311(1 1 3C(j ) 


which is defined such that the and 


^lC(j) vectors are coincident, the I3CQ) vector lies along the (3) 
or outer pivot, and the l2c(j) vector completes the orthogonal coordi- 
nate system. The base of the CMS or "Bq^" space is described by 
the ^2B(j)* an(1 ^3B(j) coor dinate system. With these spaces 

defined, the zero position of the CMG is defined when the vectors 


1 lA(j ) * 1 2A(j) ’ and 1 


are coincident with the vectors lic(j) 


l 2C (j)» and atl( * coincident with the vectors 1 2B(j) > 

andT 3B(j)' 

The (1) pivot angle ^i(j) is defined positively when the inner 
gimbal is rotated in a positive direction about the l^^jj vector with 
respect to the outer gimbal as shown in Figure III-3. Similarly, the 
(3) pivot angle is positive when the outer gimbal is rotated in 

a positive direction about the vector with respect to the base 

as shown in Figure II 1-4. 


Using the definitions of the coordinate systems and gimbal angles , 
the transformation matrices from " A (j>" t0 S P 3Ce ’ " A (3 ) " t0 

"Cqj" space, and "Cq)" to ”B^j" space can be obtained. These 
transformation matrices are 
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R ic(j) 


1 0 

0 " 


H lA(j) 

H 2C(j) 


° C6 l(j) 

" S6 Uj> 


H 2A(j) 

R 3C(j) 


0 S5 1U) 

c6 iai 


-3A(j2 


Cm- 3) 


H lB(j) 


C6 3(j) " S6 3(j) ° 


H ic(j) 

H 2B(j) 


S6 3(j) C6 3(j) ° 


H 2C(j) 

>(j] 


0 0 1 


R 3C(j) 


(m- 4) 


H lB(j) 


ca 3(j) _S6 3(j) C6 ia) S6 3a) S4 uj) 


H lA(j) 

H 2B(j) 

= 

S5 3a) c<5 3Cj) c6 ici> ’“30 ) S6 ia) 


H 2A(j) 

R 3B ( j) 


S6 lCj) ®1C3) 


“3A(j2 


(III-5) 


where 

S6 1(j ) - sln6 i(J) 

i - 1,3; j - 1,2,3 

and 

“i(J> ‘ C ° sS i(j) 

i = 1,3; j = 1,2,3 


This designation for sines and cosines of angles will be employed 
throughout this analysis . 

Referring to Figure III-2, the transformation matrices from the 
base coordinate system of each of the CMGs to the vehicle coordinate 
system are determined for the (1) CMG to be 
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■ " 

r 

1 

XV 


^y V 


T zv 


- 

L 


zv 


xv 


yv 


"ZV 


0 0 


and for the (2) CMG to be 
T xv 

V 


0 0 


and for the C3) CMG to be 


0 0 


0 0 


0 0 


0 



- — 


^IB(I) 


1 2BC1) 

- 

1 3B (!) 


r -k 


1 1B(2) 


1 2B(2) 

- 

X 3B(2) 

_ 

- — 


1 1B(3) 


1 2B(3) 


1 3B(3) 


(III-6) 


(IH-7) 


(III-8) 


Equations CIH _ 6) , (III-7) , and (III-8) can be used to determine the 
resultant momentum vector of the CMG cluster in vehicle coordinates as 


«XV 


V 

= 

H 


zv 



H 2B(1) + H 1B(2) + H 3B(3) 


H 3B(1) + H 2B(2) + H 1BC3) 


H 1BC1) + H 3B(2) + H 2B (3) 


(III-9) 
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Now, since the vector l2A(j) s P ace is defined to be along the 

spin reference axis, the transformation from "Aqj" space to "B^" 
space can be used to obtain the resultant momentum vector of the CMG 
configuration in terms of the momentum magnitude of each CMG. Assuming 
the momentum for each CMG to be Hg, the momentum vector for the j th 
CMG can be written in space as 


s (j) 


0 


_H lA(j ) " 

H o 

= 

H 2A(j) 

o 

- 


_ H 3A(j)_ 


(III-10) 


Substituting Equation CHI-10) into Equation (III- 5) gives the momentum 


l a> 


CMG in 

space 

as 

-S6 3(j) C6 l(j) 


H lB(j) 

C5 3(j) C6 l(j) 

X 

o 

li 

H 2B(j) 

S6 l(j) j 


H 3B(j) 

J 


— _ 


(in-11) 


Substituting the expressions for the CMG momentums in space from 

Equation (III-ll) into Equation (III-9) yields the expression for the 
momentum of the CMG cluster in vehicle coordinates as 


"Hxv" 


C6 3(1) C6 1C1) ■ S<5 3(2) C6 1(2) + S6 1(3) " 

H 

yv 

= 

S5 1(1) + C5 3(2) C6 1(2) _S<5 3(3) C6 1(3) 

B zv 


- S6 3C1) C5 H1) +S5 l(2) + C6 3(3) C6 1(3) 


(III-12) 



w+* 
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Equation (III-l) can be expressed in vehicle coordinates so that the 
torque applied to the vehicle becomes 


T _ 



+ fixH = H+ QxH 


vehicle space 


(III-13) 


where H is as defined in Equation (111-12) and fi is the total angular 
velocity of the spacecraft 


$ = 


X i 


(III-14) 


L- -J 

indicates the time derivative of H with respect to vehicle space. 
Substituting the expressions for ft given by Equation (III-12) and per- 
forming the indicated differentiation and cross-product operations yield 


+ 

T = D6 + Eft 


(III- 15) 


where 


" S6 3a) c6 i(i) 

~ C6 3(1) S5 1(1) 

- C6 3(2) C6 1 ( 2 ) 

0 

C5 1(1) 

- S(S 3(2) c6 1(2) 

_C5 3(1) C5 1(1) 

SS 3Ci) M iCD 

0 

S5 3(2) Sfi l(2) 

0 

C *l(3) 

~ c5 3(2) S5 1(2) 

-C5 3(3) C5 1(3) 

S5 3(3) S,J 1(3) 

C5 1(2) 

_SiS 3(3) C5 1(3) 

" C5 3(3) S6 1(3) 


(III-16) 
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and 


E = 


t = 


0 ** 

- R z 0 

«y 

fi 3Cl) 

hd) 

» 

5 3(2) 

$ 1 ( 2 ) 

• | 

*3C3) 

*1(3) 



(III-17) 


(IXI-18) 


B. Some Basic Approaches to CMG Steering Laws. 

Basic CMG steering laws are concerned with commanding gimbal rates 
that will cause the moment exerted on the vehicle by the CMGs to re- 
produce the commanded torque as nearly as possible. The general form 
of vehicle attitude control using CMGs is shown in Figure III-5. The 
analysis will consider the CMGs to be free of any gimbal stops. The 
steering laws considered in this section assume that a vehicle control 
law to generate the commanded torque is available. Simulation results 
are presented in Chapter V for the steering laws presented in this 
section. 

Since the D matrix is not a square matrix, it does not have an 
inverse and Equation (111-15) cannot be solved for exactly. The 
problem is to take any three-dimensional torque vector command T c and 
develop a six-dimensional gimbal rate vector command for the CMG 
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torque motors. As pointed out by Ross in f3j , one obvious solution 
to Equation (III-15) is 


- D@CT n - 


Efi) 


CIII-19) 


where D© represents a matrix which satisfies the following equation 


DD 




(IXI-20) 


where I is the identity matrix. There are an infinite number of D® 
matrices which satisfy this relationship. For example, partitioning 
the D matrix so that it can be expressed as 


D = {D x j D 2 J 

suggests that a possible solution to Equation (111-20) is 



( 111 - 21 ) 


(111-22) 


This solution is valid for those orientations where the D]_ and D 2 matrices 
are nonsingular. 

The above concept can be extended to the case where the solution 
to Equation (III-15) is 

i _ T -► -*■ 

d c = D ( T c - EB). (111-23) 

where D~ r is the pseudoinverse of D. D -1 can be chosen to be the pseudo- 
inverse that gives the solution for the commanded gimbal rates with 




Figure III-5. Vehicle attitude control with control 


moment gyros. 
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the minimum norm if the solution is not unique 14] . This can be 
thought of as the fi vector that has the minimum length of all possible 
solutions. This technique does not depend on any matrices being 
non-singular. 

It is desirable to design a control scheme that would minimize 
the amount of control effort required to achieve a given quality of 
control. One approach [5] is to find the law which minimizes a per- 
formance index such as 


J ” l 

i-1,3 

j-1,2,3 


S i(j) q i 


k=x,y,z 


The q^s relatively weight the control effort desired from the inner 
and outer torque motors with respect to the quality of control. In 
order to obtain an algebraic control law In terms of system parameters , 
the performance index is presented in terms of summations. This 
evaluates the system on a point-by-point basis in time which indicates 
that this law actually places a stronger constraint on the system 
than the integral performance measure which evaluates the system per- 
formance on an average basis. Expressing Equation (III-24) in vector 
form gives 

♦ a 

J = 6 T Q 5 + (T c - T) x (T c - T) (III-25) 

Substituting Equation (III-15) into (III-25) yields 
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J = £ T Q d + (x - - E^) T (T n - D$ - 


m) 


* • 


->t A -> -*T ->■ ->T ^ -KT ■* Vr T ^"T T ’*£ 

= 6 1 Q £ + T T - T D<$ - T Eft - a D T + J D 1 D5 

' r f p r> 


->-T T* -* -^r X “► -HT T ■+■ ->T T "► 

+ « D Efl - fi E T + fi E D« + fl E 1 Ef2 

c 


(III-26) 


In order to minimize J, the partial derivative 9J/9o^ must be equated 
to the zero vector. Simplifying this expression yields 


3J 

1T = 


Q t - D T (T - Efi) + D T D t 


2 = 0 


(III-27) 


Solving Equation (III-27) for 8 gives 

* c = IQ + D T D] -1 D T (T c - eS) (111-28) 

A 

Q is the positive definite weighting matrix determining the relative 

control effort of the inner and outer gimbal torque motors. Defining 

'■ 1 

Q as K_g ^ , the steering law gain, and referring to Equation (III-24) , 

it is observed that increasing K gL will increase the quality of control 

at the expense of higher gimbal rates. Figure III-6 shows a scheme 

T 

to realize this CMG control law with D referred to as the transpose 
Steering Law " I3J. 

The above steering laws have not considered the problem of gimbal 
stops or antiparallel orientations. While the laws can be modified to 
protect against these conditions, this is not considered in this 
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investigation. The basic steering laws developed tn this chapter 
from the CMG equations of motion are compared to a total optimal CMG 
control law (presented in Chapter IV) in Chapter V. 



Figure III-6. A scheme to implement the transpose steering law. 




IV. A TOTAL OPTIMAL CMG ATTITUDE CONTROL LAW 


The steering laws presented in Chapter III assumed that a controller 
provided a commanded torque as input to the CMG steering law. In this 
chapter a control law Is developed that commands the CMG gimbal rates 
based on the input of the attitude error angles and the error rates. 

The variables of the non-linear, multi-input, multi-output CMG system 
are treated in a special way so that the tools of linear system theory 
can be applied to this problem. The feedback control policy is determined 
using the state variable approach and optimal control theory. The 
development of this control policy utilizes Kalman’s work showing that 
the optimal control for a non-linear system can be given by the solution 
of the optimization problem for a set of system equations linearized 
about the current operating point [6]. 

The system presented in this chapter is a special case of the one 
developed by Skelton in [7] using the concepts described above. The 
system presented here differs in that the gyros are considered to be 
free of gimbal stops and the problem of antiparallel orientations is not 
considered. The result of applying this approach is a closed loop law 
that generates the six gimbal rate commands as a linear combination of 
the vehicle rates and vehicle attitudes. The gains multiplying these 
variables are updated as a function of the changes in the gimbal angles. 

The states of the system are chosen to be the attitude errors, 
vehicle rate errors, and gimbal angles so that 
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(IV-1) 


(IV-2) 


(IV-3) 


The gimbal rates are chosen as the control vector so that 



(IV-4) 
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For the SIXPAC CMG system, the equation for the vehicle can be written as 

9 • 

t = J* 1 [dI + eI - 8 x J • + T 1 (IV- 5) 

where T gx = the external torque applied to the vehicle 
J = the inertia dyadic of the vehicle 
and D, E, t, and ti are as described in Chapter III. 

When the primary external torque is due to gravity gradient, it 
is shown in [8] that this torque can be expressed as 

T ex - T n (t) + G(t) t (IV-6) 


where T n (t) is a time dependent gravity gradient torque evaluated at a 
nominal attitude and the elements of the G matrix are the first order 
partial derivatives of the torque with respect to 0 and evaluated at 
the nominal 0 . 

The components of the body rates can be related to the Euler 
angle rates by manipulating the transformation from inertial space to 
vehicle space. If the Euler sequence of rotations is 1, 2, 3, then the 
vehicle reference frame (v) may be related to the inertial reference 
frame (z) as 

V - te - [ 3 eJ [ 2 0] [H] z (IV-7) 


where 



0 

cos0-£ 

-sin6^ 


0 

sin0^ 

cos0^ 


(IV~8a) 
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2 


e = 



cos8 2 

0 

sin @2 

cos 0^ 
-sin 83 
0 


0 -sin0 2 

1 0 

0 cos02 

sin0 3 0 

cos 83 0 

0 1 


(IV-8b) 


(IV-8c) 


Therefore, using Equations (IV-8) i n Equation (IV- 7 ) 
0 are determined to be 

9 11 = CO30 2 cos9 3 

= cosO^sinQ^ + sin9^sin8 2 cos0^ 

®13 - sinQ^sinO^ - cosQ-^sinQ^cosQ^ 

0 2l = -cos 9 2 sin 6 3 

0^2 ~ cos0^cos6^ - sin0^sin02Sin02 
®23 = sin0^cos0 3 + cos 6^sin02sin02 

®31 = sin0 2 
9^2 = -sin©^cos02 
833 = cos0^cos0 2 
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1 2 3 

It should be noted that the 0, 0, 0, and 0 matrices are unitary 

matrices which means that their inyerse is equal to their transpose. 
This property is used below to get an expression for the vehicle 
angular velocity. The angular velocity of the vehicle can be ex- 
pressed in terms of the Euler rates as 


fl = 9^ + $2^2' + 0 3 v 3 


= lz v V 2 % v 3 J 


r«n 

10, 


0. 


(IV-10) 


where v' is an intermediate coordinate system. Equation (IV-7) can 

i 

alternately be expressed as 


1 

V 1 


V 


V 


6 13 

v 2 

= Z 1 

0 21 

+ z 2 

0 22 

+ z 3 

0 23 

v 3 


0 3 1 

- 


9 32 


0 33 


The intermediate coordinate system v ' is 


-* 

v" 


1 I T+ 

l Q]z = [ e]e v 


(IV-ll) 


(IV- 12) 


Then solving for v in terms of v* in Equation (IV-12) yields 

,1 x.->- 

v = 0[ A e ]v' (IV-13) 

Using Equations (IV-8a) and (IV-9) and performing the indicated matrix 
multiplications in Equation (IV-13) yields 
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V 


8 n 


sin03 


v 2 

v 3 

\ 

tH 

> 

II 

0 21 

e 31 

+ v 

cos0 3 

0 

* * 

+ Y 


0,oCos0, - 0.,„sin0 
13 l 12 i 

Ojjcose^ - 022 s l n 6^ 
03^cos02 - 032 s l n ®2 


(IV-14) 


Finally, the vehicle coordinate system can also be expressed as 


V 1 


1 


r 

0 


’0' 

v 2 

= v x 

0 

+ v 2 

1 

+ V 3 

0 

v 3 


0 


0 


1 

. J 


(IV-15) 


Therefore, from Equations (IV-10) , (IV-11) , (IV-14), and (IV-15) the 
expression for the vehicle rates becomes 


Qi 




ft. 


cosSgCose^ sin0- 
-cos0 2 sin0 3 cos0- 
sin0 o 0 



— — 


8 1 


8 2 

J 

_ 8 3 


(IV-16) 


Solving Equation (IV-16) for the Euler angle rates yields 


10- 


cos0 3 /cos ©2 -sin 63/ cos ©2 0 

sin0 3 cos0 3 0 

-tan 02 cos 03 tan9 2 sin0 3 1 


fti 


ft-: 


(IV-17) 


Equations (IV-5) and (IV-17) describe the dynamics of the system to 
be controlled. These equations are expressed in terms of the state 
variables defined in Equations (IV-1) through (IV-4) as 
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(IV-18) 


where W is as given in Equation (IV-17) with 83 and 63 replaced by 
X 2 and X 3 respectively, and 



Equations (IV-18) through (IV-20) are of the general form 

x = f(x,u,t) (IV-21) 


Expanding Equation (IV-21) about the nominal operating point, 

» 

-y 

X nom ^nto t ^ ie Taylor Series for x yields 
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x « f(x 


*U 


-> 

df 


nom nom 


> t,+ t 


■> 

^om 
*+■ -* 

u=u nom 


(x - x ) + 

nom aS 


& 


(u - u ) 
- -v nom 

^om 

4- + 

US5U nom 


+ g(x - X , u 
nom 


u 


nom* 




(IV-22) 


where g contains the higher order terms. 


Defining 


y = x - ^om 


and 


V = u - u 


nom 


(IV-23) 


the general form of the linearized equations is 


y = Ay + Bv 


(IV-24) 


where 


A. . 
ij 


9f i 


and B 


ij 


= l£i 

9u. 


x=x 

nom 

u=u nom 


x=x 

^nom 

u_u nom 


The steady state gimbal rate, is obtained as the equilibrium 

V -+ 

solution of Equation (IV-5) , (i.e. Q = ft = 0) , as 

%,„« ' + G °> (IV - 25 > 

where D * is the pseudoinverse of the nonsquare D matrix and Equation 
(IV-6) has been substituted into Equation (IV-5). %om cons idered 
a constant in Equation (IV-25) since the time constants of the gravity 
gradient torques are large compared to the time constants of the 
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stabilized control system. 

It is desirable to keep the average attitude and rate errors 
small as well as prevent excessive peak errors in attitude and rate. 
Therefore, the performance measure for the optimal system penalizes 
position and rate terms to avoid excessively large control efforts. 
Since gimbal stops are not considered here, there is no reason to 
penalize the gimbal angles. The performance measure is 



l-KT •> -*t -*• -+T -*■ 

1 0‘Le + ft Jfl + 6 R6 , 
e e e ef 


dt 


(IV-26) 


where 

9 e 

->■ -y -► 4 

= 0- 6 , = 6 - 0 

command nom 

and 

M 

K 

* * 

= 6 — 6 

nom 


Since ^om =0, Equation (IV-26) can be written in terms of the defined 
state variables as 



Ptu) 

J 0 

j-vT ->■ , 
W Qy + 

+T -J . 

v RvMt 

(IV-27) 

where 


L 3x3 

0, Q 

3x3 

0, . 
3x3 



Q - 

°3x3 

J 3x3 

°3x3 




_°3x3 

°3x3 

°3x3_ 




and 



The term represents the vehicle kinetic energy and minimiz- 

ing the integral of the vehicle kinetic energy conserves vehicle energy 
expenditure. L and R elements are generally chosen smaller than 
J because of this aspect of energy conservation, 

Kalman shows in [6J that a system described by Equation (IV-24) 
minimizes a functional of the form of Equation (IV-27) (where Q is any 
positive semi-definite matrix and R is any positive definite matrix) 
subject to the constraints 


y (0) “ y Q 
y (») = unspecified 


when 


V 


-1 T ■+ 
HR X B Ky 


(IV-28) 


where K is the symmetric, positive definite gain matrix that is the 
solution to the matrix Riccati equation 


K = -KA - A T K - Q + kbr _1 b t k 


(IV-29) 
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Equations (IV-28) and (IV-29) give the optimal solution if the system 
given by Equation (IV-24) is completely controllable. If the system 
is not completely controllable, the solution is optimal only when 
the uncontrollable states are not penalized. 

The solution to the linear regulator problem is the solution to 
the non-linear plant if the linearized model is updated frequently. 

The linearized model must be updated as a function of time, gimbal 
angles, and vehicle attitude commands. The gimbal angles will be the 
fastest to change and usually can provide the criterion for updating 
the linearized model. Equations (IV-28) and (IV-29) give a sub- 
optimal solution as a result of the frequent updating of the linearized 
model, but if the updating period is long compared to the system time 
constants, then this solution approaches the optimal. An implementation 
of this sub-optimal control system is shown in Figure IV-1. 

A. Sensitivity Analysis of the Riccati Equation 

A difficulty in implementing the optimal control system given by 
Equations (IV-28) and (IV-29) centers about finding the steady state 
solution to the matrix Riccati equation given by Equation (IV-29) . 
Problems encountered are excessive computation time for the desired 
accuracy of the Riccati equation solution. Two approaches are consider- 
ed to decrease the degree of these difficulties. First a sensitivity 
analysis of the gains with respect to the gimbal angles is performed 
to determine if some gains may be considered constant. Those constant 
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gains would not have to be computed at each linearization update, thus 

saving on computation time. In Section B different solution techniques 

are evaluated with respect to computation time and accuracy to determine 

which solution method is most acceptable with respect to these requirements. 

As is pointed out by Skelton in [7], it is sufficient to solve the 

sixth order Riccati equation when the gimbal angles are not penalized. 

The reduced sixth order system is 


- X, „ 


“0 

w " 


X- _"1 


'0 " 

1-3 

- 




1-3 

+ 


- x 4-6- 



A 2- 


- X 4-6- 


- B l- 


(IV- 30) 


where W is given in Equation (IV-17) 


A„ = J 


-1 


0 -h„ h 


h 3 0 


L- h 2 h l 


2 

-h. 


B. = 


1 x=x 


nom 


and J is the inertia dyadic of the vehicle. Since the gimbal angles are 
not penalized in the performance measure, only the sixth order case 
will be considered here. 

Rewriting Equation (IV-29) in the steady state and showing explicit 
dependence on a parameter a yields 


A T (a)K(a) + K(a) - K(a)S(a)K(ct) + Q = 0 


(IV-31) 
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where S(o) = B(a)R ^B T (a) 

Taking the partial of Equation (IV-31) with respect to a gives 


|^ k + a t|k + |k a+k |A.|k sk _ k |S k _ ks |k 
da 3a da da da da da 


(IV-32) 


Combining like terms simplifies Equation (IV-32) to 

^ (A - SK) + (A - SK) T H = -Kn - n T K + K6K (IV-33) 

where n = dA/da 

6 = dS/da 

and with all matrices evaluated at a - a^. The form of Equation (IV-33) 
is 

FG + G X F + H = 0 (IV-34) 

where F = dK/da 

G = A - SK 

and H = Kn + q T K - K6K 

Kleinman presents an algorithm in [ 9] that is used to iteratively solve 
equations of the same form as Equation (IV-34) numerically. The 
steady state Riccati solution for several initial gimibal angle configura- 
tions (K(ag)) is computed here; then Equation (IV-34) is iteratively 
solved for dK/da using Kleinman 1 s algorithm. The percentage variation 
in K for a Aa change in a gimbal angle is approximated as 
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q K _ dK Aa ( T v— IS - ) 

s a ^5 |T UV 

where S a is the sensitivity of K with respect to parameter a 

a is a gimbal angle, a- 6, 6 , ...6 

12 6 

and Aa is the incremental change in the gimbal angle a. 

Performing these operations for gimbal angle orientations that are 
combinations of 0° and 30° gimbal angles is sufficient to indicate 
that no gains are insensitive enough to gimbal angle variations to 
consider them constant. 

The diagonal elements of the steady-state solution to the Riccati 
equation are generally much less sensitive to gimbal angle variations 
than are the off-diagonal terms. Table IV-1 gives the range of the 
values of the off-diagonal terms of the gain matrix for the various 
gimbal angle configurations that are considered. For these same cases , 
Table IV-2 gives the range of the sensitivities of these off-diagonal 
terms to 5° incremental changes in the 6^ gimbal angle. Results are 
presented only for the 5^ gimbal angle sensitivity because these re- 
sults are typical of those obtained for the sensitivities of these 
elements to the other five gimbal angles. Since Table IV-2 indicates 
the extreme sensitivity of the off-diagonal terms to gimbal angle 
variations, only the diagonal terms can possibly be insensitive enough 
to gimbal angle variations to be considered constant. 

The range of the values of the Riccati gain matrix diagonal 
terms are given in Table IV-3 for the different gimbal angle configura- 
tions analyzed. Table IV-4 describes the range of the maximum 
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Table IV-1. Range of Riccati Gain Matrix Off-Diagonal 
Elements For Several Combinations of 0® 
and 30° Gimbal Angles. 


Gain Matrix Element 

Range of Matrix Element Values 


Max 

Min 

K 12 

1.75 x 10 4 

-4.32 x 10 4 

K 13 

5.27 x 10 3 

-4.30 x 10 4 

K 14 

2.51 x 10 5 

1.55 x 10 5 

K 15 

6.09 x 10 4 

-1.51 x 10 5 

K 16 

1.91 x 10 4 

-1.49 x 10 5 

K 23 

4.34 x 10 4 

-8.03 x 10 4 

K 24 

6.14 x 10 4 

-1.49 x 10 5 

K 25 

1.60 x 10^ 

1.28 x 10 6 

K 26 

1.55 x 10 5 

-2.86 x 10 5 

K 34 

1.83 x 10 4 

-1.49 x 10 5 

K 35 

1.56 x 10 5 

-2.85 x 10 5 

K 36 

1.84 x 10 6 

1.27 x 10 6 

K 45 

2.52 x 10 5 

-6.08 x 10 5 

K 46 

7.84 x 10 4 

-6.06 x 10 5 

K 56 

6.67 x 10 5 

-1.20 x 10 6 
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Table IV-2. Range of Sensitivity of Riccati Gain Matrix Off- 

Diagonal Elements to Variations in S. Gimbal Angle 
For Several Combinations of 0® and 30° Gimbal Angles. 
5° Variation in 6^ is Used. 


Gain Matrix 
Element 

Range of Sensitivity of Matrix Element Values 
(Absolute Value) 

, Max 

Min 

K 12 

7.77 x 10 5 

4.42 x 10 -3 

K 13 

2.43 x 10 6 

3.05 x 10" 1 

K 14 

9.81 

7.74 x 10 -4 

K 15 

2.79 x 10 4 

2.11 x 10 -2 

K 16 

7.36 x 10 

2.73 x 10 _1 

K 23 

4.33 x 10 4 

7.83 x 10~ 3 

K 24 

2.39 x 10 3 

1.23 x 10 -2 

K 25 

1.53 

2.09 x 10 -6 

K 26 

1.35 x 10 4 

1.54 x 10 -2 

K 34 

9.55 x 10 3 

1 

1.12 x 10 -1 

K 35 

8.41 x 10 3 

5.89 x 10 4 

K 36 

2.74 

1.01 x 10 -3 

K 45 

4.04 x 10 5 

3.23 x 10 -3 

K 46 

1.22 x 10 8 

4.37 x 10 -2 

K 56 

5.79 x 10 4 

1.68 x 10 3 
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Table IV-3. Range of Riccati Gain Matrix Diagonal Elements For 
Several Combinations of 0° and 30° Gimbal Angles. 


Gain Matrix 

Range of Matrix Element Values 

Element 

Max 

Min 

K 11 

3.83 x 10 5 

3.54 x 10 5 

CM 

CM 

2.39 x 10 6 

2.30 x 10 6 

K 33 

2.43 x 10 6 

2.28 x 10 6 

K 44 

8.89 x 10 5 

5.10 x 10 5 

K 55 

5.78 x 10 6 

4.41 x 10 6 

K 66 

6.84 x 10 6 

4.38 x 10 6 


Table IV-4. Range of Maximum Value of Sensitivities (for 5° Gimbal 
Angle Increments) of Gain Matrix Diagonal Elements to 
the Six Gimbal Angles. 


Gain Matrix 
Element 

Range of Sensitivity of Matrix Element Values 
(Absolute Value) 


Max 

Min 

K 11 

1.30 

2.93 x 10‘ 2 

K 22 

2.93 

1.79 x 10" 1 

K 33 

1.47 

1.93 x 10” 1 

K 44 

4.66 

3.94 x 10* 1 

K 55 

5.86 

1.25 

K 66 

4.81 

1.27 
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sensitivities of these gain elements with respect to variations in 
the six gimbal angles. Table IV-4 illustrates that the diagonal terms 
are much less sensitive to the gimbal angles than the off-diagonal 
terms; however, it also indicates that no elements of the Riccati gain 
matrix are insensitive enough to the gimbal angles to be considered 
constant. As a consequence, all elements of the gain matrix must be 
computed at each linearization update and no computational time can 
be saved by considering some elements to be constant. 

These results were obtained by using the same system as Skelton 
in [ 7J. Therefore, the procedure employed in this section should be 
employed for the system under consideration. If an element of the 
gain matrix is determined to be insensitive to given parameter varia- 
tions within a specified degree, then that element may be considered 
constant thus reducing required confutational time. 

B. Numerical Solution Techniques for the Riccati Equation 

A conventional Runge-Kutta integration scheme is used in [ 7] to 
solve Equation IV-29 for K in the steady state. Some other numerical 
techniques are considered here and compared to determine if another 
algorithm might give the Riccati gain matrix with less computational 
time and accuracy comparable to a Runge-Kutta method. 

Two other integration algorithms, Euler and Modified Euler, are 
evaluated along with the Runge-Kutta scheme. These should have smaller 
computational times at a cost to the accuracy of the steady state solution. 
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Russell presents a mathematical technique in [ 11] that accelerates 
the convergance of a transient computer solution to steady state con- 
ditions thus substantially reducing execution time. This procedure is 
applicable when 

(1) the boundary conditions and internal constraints are 
either constant or cyclic 

(2) the steady state results are independent of the initial 
conditions , and 

(3) only the steady state solution is desired. 

There are two accelerations techniques — exponential extrapolation 
and constant proportionality adjustment. Exponential extrapolation 
involves fitting an exponential curve through three points in the 
transient solution and predicting the steady state value. The general 
form of the fitting curve is given as 


Y(t) = + C 2 e“ C 3 t 


(IV- 36) 


Assuming the values of the three points along the transient solution 
to be Y^, Y 2 , and Y^, then the steady state value for Y(t) in Equation 
(IV-36) is given as 


Y (t 


8 ) = 


- y 2 

? 1 + *3 ' 2I 2 


(IV-37) 


The solution adjustment is obtained from Equation (IV-37) as 

2 

Y Y - Y 

AY = — L-3 2 _ Y 

Y 1 + Y 3 - 2Y 2 1 


(IV- 38) 
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The constant proportionality adjustment is 


AY = SKIP 



At 


(IV- 39) 


where SKIP is a constant multiplication factor input by the user. 

This method is usually less accurate than exponential extrapolation, 
but it usually requires less computation time as well as requiring 
less storage space. 

Since the possibility of solution instability is enhanced with 
the accelerated integration techniques, the sign of the adjustment to 
the gain values can be checked. If the sign changes, then the in- 
cremental adjustment is modified by a constant reduction factor to 
prevent instability or at least to retard it. If the sign change is 
indicative of a simple overshoot rather than instability, the modifica- 
tion to the incremental adjustment should speed up the solution computa- 
tionally. 

In [ 9J and [ 12J Kleinman presents an iterative algebraic technique 
to solve for the linear regulator Riccati gain matrix. If V^, k = 0, 

1 , ... is the unique positive definite solution of the linear algebraic 
equation 

v k + \\ + Q + RL k - 0 (IV-40) 

whe re re curs i vely , 

\ = R _1 B T V kl k = 1, 2, ... (IV-41) 

\ = A - BL k 


(IV-42) 
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and where Lq is chosen such that the matrix Aq = A - BLq has eigenvalues 
with negative real parts, then 


K£\ +1 <\£... k=0, 1 


(IV-43) 


and 



CO 



(IV-44) 


where K is the steady state solution to Equation (IV-29). Kleinman 
presents the computational algorithms to implement this procedure in [9]. 

A program for the Linear Quadratic Loss (LQL) problem ( 13] uses 
Potter's algebraic method [14] to obtain the steady state solution to 
the Riccati equation. Potter’s method involves finding the eigenvectors 
(or pseudo eigenvectors) corresponding to eigenvalues with negative 
real parts of the 2n x 2n Hamiltonian matrix. 


-1 T 

A -BR B 



(IV-45) 


Spectral factorization of the Hamiltonian matrix is used to obtain these 
eigenvalues. The stable eigenvectors are then used to form a 2n x n 
matrix whose columns are the real eigenvectors. If the eigenvectors 
are in complex conjugate pairs, two vectors made up of the real and 
imaginary parts of one of the complex eigenvectors are used instead 
of the complex entries. If this 2n x n matrix so formed is 


D 

E 


(IV-46) 
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then the steady state solution to the Riccati equation is determined 
to be 

K = ED _1 (IV-47) 

Using the same system parameters as in [ 7] , solutions to the 
Riccati equation are obtained using the previous algorithms. Since 
the algorithms would be programmed on a particular flight control com- 
puter, the actual computation times for these procedures would be 
meaningless . Therefore, the tabular results presented in the remainder 
of this chapter are normalized. The minimum value of the parameter for 
all cases considered is assigned the value zero while the maximum 
value is assigned the value 100. All intermediate values are assigned 
a value between 0 and 100 in the following linear fashion, 

— — EsiE 100 (IV-48) 

Pmax Pmin 

where p n is the normalized parameter value, 

p^ is the actual parameter value, 

Pmin is tiie parameter value, 

and p is the maximum parameter value. 

L max 

Since the steady state solution to the Riccati equation is desired, 
the derivatives of each of the gains at the computed solution point 
gives a measure of the solution accuracy. For the results included in 
the remainder of this chapter, the accuracy comparison values are 
obtained by summing the absolute values of all the elements of the 
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gain derivative matrix. These values are then normalized according 
to Equation (IV-48). 

Euler, Modified Euler, and Runge-Kutta integration procedures 
are examined in this analysis. Additionally, these solution techniques 
are examined with exponential acceleration and proportional accelera- 
tion, both with and without the instability retardation routine. 
Finally, the Kleinman and Potter algebraic routines are compared to 
the integration routines. Since a conventional Runge-Kutta integra- 
tion algorithm is used in [ 7] , this case is included in several of the 
following tables for the purpose of comparison. 

The integration routines are very sensitive to the integration 
time increment. Table IV-5 illustrates the expected decrease in com- 
putational time for an increase in the integration time increment for 
the conventional Runge-Kutta integration solution. 

For those cases including instability suppression, the value of 
the constant factor to modify the solution adjustment is chosen to be 
.1. This value represents a good choice, but not necessarily a best 
choice. Solutions to the Riccati equation are obtained using other 
values, but there is no general trend. In some instances .1 gives the 
better result with respect to computational time and accuracy, while 
in other instances .025 and .25 give better results for the three 
values. Table IV-6 compares the three integration methods to each 
other with and without instability suppression. The instability 
suppression has a more noticeable effect on the Euler and Modified 
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Table IV-5. Computational Time and Accuracy for the 
Runge-Kutta Integration as a Function of 
the Integration Time Increment (At). 


At 

Normalized 

Accuracy 

Normalized Computational 
Time 

-.0001 

* 

* 

-.005 

* 

* 

-.01 

* 

* 

-.05 

2.79 

i 

36.17 

-.1 

1.59 

17.86 

-.5 

.14 

3.19 


* indicates that the total computational time was larger 
than that allowed for this analysis. 


Table IV-6. Comparison of Integration Solution 
Methods For At = -.5. 


Solution Method 

• 

Normalized 

Accuracy 

Normalized 

Computational 

Time 

Normalized 
Weights of 
Time & Accuracy 

Runge-Kutta 

.14 

3.19 

3.33 

Runge-Kutta 

(Instability Suppression) 

.14 

3.23 

3.37 

Modified Euler 

.07 

2.91 

2.98 

Modified Euler 

(Instability Suppression) 

.17 

2.23 

2.40 

Euler 

0 

.08 

.08 

Euler 

(Instability Suppression) 

.10 

7.77 

7.87 
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Euler routines. If computational time and accuracy are weighted 
equally as in Table IV-6, then the Euler algorithm accomplishes the 
desired results most efficiently. 

The exponential acceleration of integration solutions is discarded 
for this system because the transient solution is not approximated by 
an exponential very well. As a result, the solution either goes un- 
stable or requires a prohibitive amount of computation time if the 
instability suppression is included. 

The linearly accelerated integration algorithms are sensitive to 
the value of the constant of proportionality used to accelerate the 
soltuion. If this constant is chosen too large, this procedure will 
also result in unstable solutions or prohibitive computation times if 
instability suppression is used. Table IV-7 provides a comparison of 
some solutions to Equation (IV-29) obtained by Euler integration and 
linearly accelerated Euler integration. The linearly accelerated 
routine produces slightly better results than the simple Euler integra- 
tion routine. For a particular system, a value for the constant pro- 
portion adjustment can be obtained empirically. This must be done to 
minimize the possibility of unstable solutions. 

For all solution techniques considered. Table IV-8 presents those 
algorithms requiring the least computational time, Table IV-9 presents 
those techniques giving the most accurate solutions, and Table IV-10 
presents the best solution methods for an equal weighting of time and 
accuracy. The Euler integration routine with a step-size of -.5 gives 
the best results of the algorithms employed to solve the Riccati 
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Table IV-7. Comparison of Euler Integration and 

Linearly Accelerated Euler Integration 
Solutions to the Riccati Equation. 


Algorithm 

At 

Normalized 

Accuracy 

Normalized 
Comp ut at ional 
Time 

Equal Weighting 
of Time and Accuracy 

Euler 

B 

1.01 

3.54 

4.55 

Euler 


2.98 

8.22 

11.19 

Euler 

1 

29.65 

34.58 

64.23 

Acc Euler 

B 

.08 

3.05 

3.14 

Acc Euler 

-.05 

1.93 

2.51 

4.43 

Acc Euler 

-.01 

1 

32.04 

12.41 

44.44 


Table IV-8. The Solution Methods Requiring the Least Computational Time. 


Algorithm 

At 

Normalized Computational Time 

Potter's Algebraic 

* 

0 

Euler 

-.5 

.08 

Euler (Linear Acceleration) 

-.1 

1.65 

Modified Euler 

-.5 

2.23 

Euler (Linear Acceleration) 

-.05 

2.51 

Runge-Kutta 

(for conq>arison) 

-.5 

3.19 


* An Algebraic Technique 
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Table IV-9. The Solution Methods Resulting in the Best Accuracy. 


Algorithm 

At 

Normalized Accuracy 

Euler 

-.5 

0 

Modified Euler 

-.5 

.07 

Euler (Linear Acceleration) 

-.1 

.08 

Euler (Instability Suppression) 

-.5 

.10 

Runge-Kutta 

-.5 

.14 

Kleiaman Iterative 

* 

.17 


* Algebraic Algorithm 


Table IV-10. The Solution Methods Providing the Best 
Tradeoff Between Time and Accuracy. 


Algorithm 

At 

Normalized Weighting of 
Accuracy and Con?) Time 

Euler 

-.5 

.08 

Potter's Method 

* 

.95 

Euler (Linear Acceleration) 

-.1 

2.22 

Modified Euler 

' -.5 

2.40 

Euler (Linear Acceleration) 

-.05 

3.28 


* Algebraic Algorithm 
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equation. Table IV-11 compares the two algebraic techniques with the 
Euler integration solution because the integration solutions have the 
property that the numerical solution goes unstable if the step-size is 
chosen to be too large. In contrast, larger step-sizes give better 
accuracy and computational time trade-offs. This is not true for the 
algebraic techniques so that, depending on the application, it may be 
desirable to use the less efficient algebraic algorithms. 
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Table IV-11. Algebraic Techniques Compared With 
the Euler Algorithm With At = -.5 
and the Runge-Kutta Algorithm With 
At = -.5. 









V. SIMULATION RESULTS 


Simulation results are presented in this chapter to compare the 
basic CMG steering laws presented in Chapter III with the total 
optimal CMG control law presented in Chapter IV. The system parameters 
used in these simulations are those used in [7]. 

It is stated in Chapter III that the three basic steering laws 
presented there require a controller that provides a commanded torque 
as an input to these steering laws. Since these basic steering laws 
are compared to the total optimal control law, the torque commanded 
by the controller will be determined from optimal control considerations. 
These conditions should provide a fair comparison of the system 
responses. 

The control law that gives the commanded torque for the basic 
steering laws is developed in a manner analogous to that used to 
develop the total optimal control law in Chapter IV. Equation (IV-27) 
again describes the performance measure, but the terms in the per- 
formance measure refer to different variables. The sixth order state 
vector is condos ed of the vehicle attitude errors and attitude rate 
errors while the control vector is the torque command vector. The R 
matrix is chosen as the identity matrix while the Q matrix relatively 
weights the elements of the state vector as well as relatively weighting 
the state vector with control effort expenditure. Linearizing the 
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state equations about the commanded vehicle attitude and the nominal 
vehicle rates puts this problem in the form of the linear regulator 
problem as in Chapter IV. The resulting optimal control torque is 
linear negative feedback of the position and rate errors. If a par- 
ticular response is desired, the closed-loop system matrix may be 
examined and varied by varying the state feedback to obtain the desired 
closed-loop eigenvalues since the optimal law is known to be constant 
negative feedback of the states. This approach will result in the 
desired response without being concerned with the problem of determin- 
ing the weighting matrices to give this response. 

The partitioned inverse steering law described by Equation (III-19) 
and the pseudoinverse steering law described by Equation (III-23) are 
both open-loop steering laws. (Algorithms for computing the pseudo- 
inverse are given in [15]). As a result, these steering strategies do 
not consider the CMG cross-coupling effects. Consequently, the system 
using these steering laws goes unstable as shown in Figure V-l. This 
response is typical for both of these open-loop steering laws. 

These steering laws are modified by negatively feeding back the 
actual torque on the vehicle to force the basic steering law to create 
a torque that more closely follows the commanded torque. This negative 
feedback may be implemented by using the measured gimbal angles and 
gimbal rates to compute the actual torque using 

T = D<5 + ESI 


(V-l) 
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since D and E are both functions of the gimbal angles only. Implement- 
ing this negative feedback has the advantage of decreasing any detri- 
mental cross-coupling effects as may be observed in Figures V-2 and 
V-3. Figures V-2 through V-5 present the responses of the different 
systems to initial position and rate errors in one axis only. The 
cross-coupling effects for the control laws may be compared to one 
another by examining the axes that had no initial errors in them. 

Figure V-6 compares the responses of the four systems to an 
initial position and rate error in one axis only. The other axes are 
not included in Figures V-6 through V-ll because Figures V-2 through 
V-5 have illustrated how small the cross-coupling effects are. The 
pseudoinverse steering law and the optimal steering law give almost 
the same response. This response is underdamped, but faster than the 
partitioned inverse steering law and the total optimal CMG control law. 
The partitioned inverse steering law gives an overdamped response that 
is faster than the total optimal control law but slower than the other 
two steering laws. This partitioned inverse steering law gives the 
best response for the weighting matrices chosen in these cases. The 
response of the total optimal CMG control law is similar to the 
partitioned inverse steering law, but it is slower in achieving the 
desired control. Figures V-7 through V-ll compare the responses of 
the four systems to one another for initial errors in two axes and 
three axes. Examination of these cases reveals that the same comments 
apply to them that applied to Figure V-6. 
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Figure V-4. Attitude vs. Time for an initial error in one axis using the optimal steering law. 




Figure V-5. Attitude vs. Time for an initial error in one axi9 using the total 
optimal CMG attitude control Law. 
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Figure V-6. Comparison of the four CMG control laws for an initial error in one axis* 
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Figure V-7. Comparison of the four CMG control laws for an initial error in two axes. 
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Figure V-10. Comparison of the four CMG control laws for an initial error in three axes. 
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Figure V-ll. Comparison of the four CMG control laws for an initial error in three aoces. 
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These results are not presented to imply that one of these 
systems is better than the other in an absolute sense. The responses 
of these systems may be varied by changing the relative weightings of 
some of the parameters. The partitioned inverse steering law and the 
pseudoinverse steering law will have different responses only if the 
relative weighting of vehicle position and rate errors against control 
effort is changed in the development of the vehicle control law. 

The optimal steering law response may be varied by changing the relative 
weighting between the gimbal rates and the torque error in the develop- 
ment of the optimal steering law. It may also be varied by changing 
the vehicle control law as for the other steering laws. The response 
of the total optimal CMG control law may be changed by changing the 
relative weighting of position and rate errors of the vehicle against 
the commanded gimbal rates. Keeping in mind this possibility of varying 
these responses, the main point to be made by these simulation results 
is that the basic steering laws with an optimal vehicle control law 
can control the vehicle as well as the total optimal CMG attitude 
control law. Furthermore, these basic steering laws are implemented 
without the computational problems associated with the total optimal 
control law. It is suggested that future investigations compare and 
analyze these systems in more detail. 
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